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We test the elasticity of granular aggregates using increments of shear and volume strain in a 
numerical simulation. We find that the increment in volume strain is almost reversible, but the 
increment in shear strain is not. The strength of this irreversibility increases as the average number 
of contacts per particle (the coordination number) decreases. For increments of volume strain, 
an elastic model that includes both average and fluctuating motions between contacting particles 
reproduces well the numerical results over the entire range of coordination numbers. For increments 
of shear strain, the theory and simulations agree quite well for high values of the coordination 
number. 
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Granular materials have received the attention of many 
researchers in the last decade because of unsolved prob- 
lems with direct relevance to chemistry, physics and en- 
gineering. The behavior of a granular material can range 
between that of a gas and that of a solid, depending on 
the applied loading and the regime of deformation con- 
sidered. Significant progresses have been made with the 
introduction of numerical tools (e.g. [l|) that permit the 
detailed analysis of an aggregate of particles. Although 
such simulations have provided information about inter- 
particle interactions, such as their elasticity, sliding, and 
deletion, and statistical measures of their cooperative be- 
havior, such as induced anisotropy and force chains, it is 
still unclear how to incorporate these informations in a 
predictive theoretical model. 

Attempts to do this have been made in the context of 
the effective medium theory (EMT) (e.g. HQ) in which 
the contact displacements are given by the applied aver- 
age strain. However, the predicted shear and bulk mod- 
uli are far from those measured in numerical simulations 
(e.g. [3, [Ej]). In order to improve the theoretical predic- 
tion, particle displacements are given by the sum of an 
average and the fluctuation components (e.g. [f| 0|). In 
particular, better predictions of the shear and bulk mod- 
uli have been obtained in a recent work that em ploy s pair 
fluctuations [!,[!]■ However, Agnolin and Roux |lfjj point 
out that such predictions fail when the the coordination 
number, Z , is low and close to the isostatic limit Zi SO at 
which the aggregate is statically determinate. They sug- 
gest that for low coordination number, more complicated 
models are needed that account for collective deforma- 
tions among particles. That is, the assumption of pair 
fluctuations is not sufficient to capture the response of 
poorly coordinate aggregates of particles. Here, we ad- 
dress this issue and provide an alternative interpretation 
of the failure of an elastic description at low coordination 
numbers. 



We carry out numerical simulations using a distinct 
element method and focus on the first incremental re- 
sponse of an isotropically compressed random aggregate 
that consists of identical, elastic, frictional spheres. We 
consider dense aggregates, with solid volume fractions <6 
near 0.64, that have different Z. Previous work [10Hl2j 
have considered poorly coordinated aggregates and found 
that the shear modulus is proportional to Z — Zi SO , where 
Zi SO is equal to four for a packing of frictional spheres. 
Here, we find that when Z decreases, there is an irre- 
versible behavior of the aggregate that involves local, co- 
ordinated, irreversible motions of the particles that are 
not resisted by forces. These motions result in a reduc- 
tion of the apparent stiffness of the aggregate (e.g. [H[). 
That is, the initial configuration of a poorly coordinated 
aggregate can not sustain any incremental strain unless 
a change in the geometry of the packing occurs. When 
such irreversible changes are present, simple elastic the- 
ory can not reproduce the response of the aggregate and 
the utility of the elastic moduli is questionable. 

We introduce measures of these irreversible deforma- 
tions which vary with the coordination number. For in- 
crements in volume strain, the number of irreversible mo- 
tions is so small that their effect is negligible and the 
response of the aggregate can be assumed to be elastic. 
More importantly, for shear increments, the strength of 
the irreversibility persists even at high Z and increases 
as Z decreases towards its isostatic value, indicating that 
elasticity does not describe the aggregate response. 

Our numerical simulations consider 10, 000 particles, 
each with diameter d = 0.2 mm, randomly generated 
in a periodic cubic cell. We employ material properties 
typical of glass spheres: a shear modulus fi = 29 GPa 
and a Poisson's ratio, v = 0.2. The interaction between 
particles is a non-central contact force in which the nor- 
mal component is the non-linear Hertz interaction and 
the tangential component is bilinear: an initial elastic 
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FIG. 1: Measurement of irreversibility in terms of displace- 
ments when increments of volume and shear strain are ap- 
plied. 
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FIG. 2: Measurement of irreversibility in terms of contact 
forces when increments of volume and shear strain are applied. 



displacement followed by Coulomb sliding (e.g. [a]). We 
create different initial isotropic states by varying the co- 
efficient of friction between particles during the prepa- 
ration (e.g. [3); ai l initial states have a solid volume 
fraction 4> ~ 0.64 and a confining pressure that varies 
from hOKPa to lOMPa. 

For all packings, we evaluate the response of the ag- 
gregate to homogeneous increments in volume and shear 
strain by setting the particle coefficient of friction high 
enough to prevent sliding. Because of the random po- 
sitions of the particles and their different initial contact 
stiffnesses, the subsequent particle motions are the sum 
of the homogeneous applied strain and a fluctuation that 
relaxes the particles towards a new equilibrium state. 
When this relaxation involves a rearrangement of par- 
ticles, the incremental response is irreversible. 

We attempt to characterize the particle rearrangement 
and introduce two measures of the strength of the irre- 
versibility: the first is related to contact displacement, 



and the second is associated with contact forces. We 
apply a forward increment in strain followed by an iden- 
tical backward strain in which the reference configura- 
tion should be recovered if the deformation were perfectly 
elastic. 

We define the parameter £ as the average over all con- 
tacts of the ratio of the absolute values of the contact 
displacements after and before the backward increment 
in strain. A parameter, x-i is similarly defined in terms 
of the contact forces. Both parameters would be zero for 
perfectly elastic behavior. 

We measure £ and x i n all °f f ne packings for incre- 
ments in both shear and volume strain. We apply incre- 
ments of strain with magnitudes that depend on the con- 
fining pressure; the ratio of the associated volume strain 
with the confining pressure is constant at about 10 ~ 2 . 
The results are plotted in Figs. [T]and[2] displaying ( and 
X as functions of Z. 

Figs. [T] and [2] show small variations in £ and x asso- 
ciated with increments in volume strain. The strength 
of the irreversibility is almost negligible, with slight in- 
creases as Zi SO is approached. We conclude that an ap- 
proximate elastic response for the aggregate is obtained 
for increments of volume strain, independent of the co- 
ordination number. 

For increments in shear strain, both £ and x increase 
as Zi SO is approached. We believe that these irreversible 
motions are associated with the presence of local insta- 
bility (e.g. [H|). Moreover, both £ and x are non-zero 
for high values of the coordination number, in contrast to 
what we find for increments in volume strain. The aggre- 
gate seems to experience rearrangements over the entire 
range of the coordination number, with the irriversibility 
becoming stronger as isostaticity is approached. 

The condition of isostaticity has been an object of great 
interest for many researchers (0, [111, UH, [l5l - tl7j V The 
contact forces in an aggregate are uniquely determined 
in terms of the applied loads, independent of the con- 
tact stiffness, when the aggregate is both statically and 
kinematically determinate (e.g. [18]). The condition for 
static determinacy, often referred to as Maxwells con- 
dition [l9|, insures the equality between the number of 
unknowns and the number of equilibrium equations. In 
a granular aggregate, this necessitates that Z = Zi SO , 
The condition for kinematic determinacy insures that 
there are no inextensional mechanisms in the aggregate; 
then a rigid aggregate is able to sustain any external self- 
equilibrated perturbation without changing the relative 
positions of its particle centers. 

The kinematic condition is not often taken into account 
in recent work on granular aggregates (e.g. [20]) and 
sometimes has been emphasized in a different way. For 
example, Moukarzel [131 ] , in his description of network 
rigidity, defines an isostatic system to be one in which 
the rank K of the rigidity matrix that relates the contact 
forces to the applied forces always equals to the number of 
equations, rather than one in which the simple Maxwell 
condition is satisfied. Then, when K is less than the 
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number of equations, the network is flexible; this can 
occur when Z = Z iso . 

The presence of irreversibility in all of our aggregates, 
even in the limit that Z — 4 (see Figs. [1] and [5]) indi- 
cates that the kinematic condition does not hold. Con- 
sequently, inextensional mechanisms and the associated 
soft or floppy vibrational modes are possible. This situa- 
tion may occur whatever the value Z , if particles location 
do not correspond to those of a rigid network. In par- 
ticular ,in packing characterized as isostatic by Z = Z iso 
(e.g. [ll|, [H, HH), soft modes can occur as long as the 
kinematic condition is not satisfied. However, we should 
emphasize that the initial states that we employ are con- 
structed in a much different way than those constructed 
to insure an initial, stable, elastic response (e.g. plj]). 

We now turn from irreversibility to elasticity and re- 
port results from numerical simulation for the bulk mod- 
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ulus, = — Aau/AV, and the shear modulus, G = 

3 i=i 

A<7i2/Aei2, where er^ is the average stress tensor, 
the average strain tensor, and V is the volume of the 
aggregate. When a shear increment is applied to poorly 
coordinated systems, the irreversibility increases and the 
material response deviates from elasticity. That is, an 
elastic theory can only be considered as providing an up- 
per bound on the shear modulus. 

To make predictions of the moduli, we adopt the fluc- 
tuation theory developed by [9(. This theory improves 
upon the simplest average strain models @, y] , because 
contacting particles are assumed to move with both the 
average deformation and fluctuations and because the 
statistics of the aggregate are taken in account. The the- 
ory employs force and moment equilibrium for a typical 
pair of particles to evaluate the fluctuations, and then 
uses them to determine the stress in the aggregate. How- 
ever, at low values of the coordination number, the simple 
statistical model introduced to describe the variability of 
the neighborhood of contacting pairs of particles is prob- 
ably too simple, ( e.g. [Hjj) and the initial distribution 
contact forces (e.g.[2J|) is not taken into account. 

Here we repair the second deficiency and assume that 
the distribution, w(P), of the normal component of the 
contact force, P is exponential (e.g.[4|): 



w{P) = i exp 



(1) 



where, by definition, P 



Pw (P) dP. 



From the effective medium theory (e.g.[24|), the con- 
fining pressure po can be expressed as function of P, 

po = Z^P/nd 2 , 

and the relation between the normal component of the 
contact displacement, <5, and P is 
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FIG. 3: Comparison between the numerical data and fluctu- 
ation theory for the normalized bulk modulus. 

Using © in ©, we obtain 
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where T is the Gamma function. 
So the average normal contact stiffness Kn 
lid 1 / 2 8 1 / 2 1 (1 — v) becomes 
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The average shear stiffness is Kt = 
2^7(1 -i/)/ (2 -i/). 

Taking in account the initial distribution of forces in 
this way, we obtain a new solution for the fluctuations; 
with this, the resulting expressions for the bulk modulus 
O and the shear modulus G are, respectively, 
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where a relation between the rms value of the fluctua- 
tion in the number of contacts per particle and its av- 
erage value has been adopted At the upper limit 
of validity of equations (|5I6[) . Z — 22/3, we recover the 
average strain prediction [2|, |3j . This limit results from 
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FIG. 4: Comparison between the numerical data and fluctu- 
ation theory for the normalized shear modulus. 



the modeling of the statistical distribution of particles in 
the assembly. 

Comparison between the predictions of the bulk and 
shear moduli and measurements in the numerical simu- 
lation are shown in Figs. [3] and 01 There is agreement 
at high values of the coordination number for the shear 
modulus; the slight difference may be attributed to the ir- 
reversibility. However, when the comparison is extended 



to lower values of Z, the prediction strongly deviate from 
the simulations, as also seen by others [10] ■ Here, we 
conclude that this discrepancy is due to the observed ir- 
reversible motions and that an elastic theory is not able 
to capture the behavior of the system for low values of 
Z. For the bulk modulus, the theory works quite well 
over the entire range of Z, because the irreversibility, 
measured by x an d C) can be neglected. 

In conclusion, we have measured the irreversible mo- 
tions in a granular aggregate subjected to incremental 
strains and determined their influence on the mechanical 
response of the material. We found that particles may 
experience rearrangements in their geometry and contact 
forces, even when small perturbations are applied and 
that these rearrangements are sensitive to the coordina- 
tion number. The strength of this irreversibility is neg- 
ligible for increments in volume strain, while it strongly 
increases for increments in shear strain as Z approaches 
Zi sa . The presence of irreversible motions in the aggre- 
gates indicates a deviation from elastic behavior. That is, 
an elastic theory is appropriate to describe the material 
behavior only if mechanisms do not play an important 
role in the displacements of the particles. 
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